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One of the distinctive features of the QCD phase diagram is the possible emergence of a critical 
endpoint. The critical region around the critical point and the path dependency of the critical 
exponents are investigated within effective chiral (2 + l)-flavor models with and without Polyakov 
loops. Results obtained in no-sea mean-field approximations where a divergent vacuum part in the 
fermion-loop contribution is neglected are compared to the renormalized ones. Furthermore, the 
modifications caused by the backreaction of the matter fluctuations on the pure Yang-Mills sys- 
tem are discussed. Higher-order, non-Gaussian moments of event-by-event distributions of various 
particle multiplicities are enhanced near the critical point and could serve as a probe to determine 
its location in the phase diagram. By means of a novel derivative technique higher-order general- 
ized quark-number susceptibilities are calculated and their sign structure in the phase diagram is 
analyzed. 
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I. INTRODUCTION 

Two important properties of the QCD vacuum are 
spontaneous chiral symmetry breaking and color con- 
finement. It is expected that chiral symmetry can be 
restored and that a confinemcnt/deconfinement transi- 
tion occurs at high temperature and/or baryon densities 
PQ. A deeper understanding of how these phase transi- 
tions manifest themselves both theoretically and experi- 
mentally is of utmost importance in mapping the phase 
diagram [2]. 

A crossover transition at vanishing density seems to 
be well established by recent QCD lattice simulations 
at almost physical quark masses. With a better con- 
tinuum extrapolation and improved actions the latest 
value of the chiral critical temperature T x w 154 ± 9 
MeV by the HotQCD Collaboration [3] is in agreement 
with T x w 147 - 157 MeV of the Wuppertal-Budapest 
group [4 . However, the extrapolations of lattice simu- 
lations towards finite chemical potential are much un- 
der debate. Most theoretical model studies evidence a 
genuine second-order critical endpoint (CEP) where the 
first-order chiral phase transition line at large chemical 
potential terminates [5J |B] . Many properties of the end- 
point such as, e.g., its exact location remain unknown. 
QCD lattice simulations cannot directly access the rele- 
vant region in the phase diagram but extensions towards 
finite chemical potential provide some limitations and can 
rule out the existence of a CEP for small fi/T ratios 
[7j. This observation agrees with recent first-principle 
QCD studies performed with functional renormalization 
group techniques and Dyson-Schwinger equations [BlU0|. 
Moreover, the predicted chiral first-order transition at 
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small temperatures is model dependent. For example, 
the inclusion of the 't Hooft anomaly in chiral mod- 
els may change the first-order transition to a smooth 
crossover [TT] . 

This underlines the importance of the planned facil- 
ities such as the Compressed Baryonic Matter (CBM) 
experiment at the Facility for Antiproton and Ion Re- 
search (FAIR) and the Nuclotron-Based Ion Collider Fa- 
cility (NICA) at the Joint Institute for Nuclear Research 
(JINR) and of the recently started "beam energy scan" at 
the Relativistic Heavy Ion Collider (RHIC, Brookhaven 
National Laboratory) [TJ] which are dedicated to search 
for a critical point. For an experimental overview see 
e.g. [IS]. The knowledge of the characteristic signatures 
of the conjectured CEP is inevitable for the experiments 
[HI [T5] . Possible experimental signatures which are 
based on the singular behavior of thermodynamic func- 
tions near a critical point were already suggested a decade 
ago (see, e.g., [16]). They are related to temperature and 
chemical potential fluctuations in event-by-event fluctu- 
ations of various particle multiplicities [17] . By scanning 
the center of mass energy and thus the baryochemical 
potential an increase and then a decrease in the number 
fluctuations of, e.g., pions and protons should be seen as 
one crosses the critical point. The nonmonotonic behav- 
ior in the number fluctuations in the vicinity of a critical 
point might serve as a probe to determine its location 
in the phase diagram assuming that the signals are not 
washed out by the expansion of the colliding system (18) . 

In a realistic heavy-ion collision the correlation length 
is cut off by critical slowing down and finite volume ef- 
fects. In addition, the system has only a finite time to 
build up correlations. Estimates of the largest correla- 
tion length as the collision cools past the critical point 
are in the range of a few (2 — 3) fm and only a factor 3 
larger than the natural scale (0.5 — 1) fm far away from 
the critical point [Tj5]. Hence, if the correlation length 
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increases only by a few percent near a critical point more 
sensible quantities are needed for the analysis of freeze- 
out and critical conditions in heavy- ion collisions [2"01l21| . 
Possible candidates are higher-order cumulants or ratios 
of higher-order generalized susceptibilities which depend 
on higher powers of the correlation length £. An example 
is the fourth-order susceptibility \4 which scales like £ 7 
near the critical point; hence, an enhancement by a fac- 
tor 2 is expected for this quantity. The sign of the fourth 
moment should also change as the CEP is approached 
from the crossover side as pointed out in Ref. [2"2"] . 

In this work we investigate the sign structure not only 
of the fourth moment (kurtosis) but also of higher-order 
moments in the (T, fi) phase diagram. We discuss the 
generic behavior of these higher moments at nonvanish- 
ing chemical potential within (2 + l)-flavor quark-meson 
(QM) and Polyakov-quark-meson (PQM) models which 
we briefly recapitulate in Sec. [TTJ The no-sea mean-field 
findings are compared to the results obtained in the 
renormalized models. As mentioned in [2"3l |2~4] the criti- 
cal region around the critical point in the phase diagram 
is not pointlike but has a richer structure. In Sec. |III| 
we investigate the size of the critical region and the path 
dependency of the critical exponents towards the crit- 
ical point. The influence of the backreaction of matter 
fluctuations on the pure Yang-Mills system on the critical 
region is analyzed. In Sec. |IV| different ratios of moments 
of the baryon number fluctuations to high orders are cal- 
culated by means of a novel derivative technique |25j . We 
conclude and point out some experimental consequences 
from these findings in Sec. W\ 



II. THREE-FLAVOR MODELS 

Chiral symmetry and its spontaneous breaking for 
three flavor is very well described by renormalizable 
quark-meson type models which serve as an effective re- 
alization of low-energy strongly- interacting matter [2"r?] . 
These models can be augmented with the Polyakov loop 
(see, e.g., (27J [28]), yielding the PQM models which 
mimic in addition to the chiral symmetry breaking cer- 
tain aspects of a statistical confinement. 

The Lagrangian of the three-flavor (Nf = 3) PQM 
model [29] is assembled from three parts 

£pqm = q{ip+h(t>5 + lo^f)q + L m -U(®,$) , (1) 

where the first term describes the interaction between the 
N c = 3 color quark fields q = (u, d, s) and the mesons 
<ps = T a (<7 a + 175 7r a ) with a flavor-blind Yukawa cou- 
pling h. By T a the nine generators of the E/(3) sym- 
metry are denoted and a a (w a ) labels the (pseudo)scalar 
meson nonets. For each flavor / G {u, d, s} an inde- 
pendent quark chemical potential /i/ is introduced but 
in this paper we will consider symmetric quark matter 
with one uniform quark chemical potential fi = /j.b/3- 
The quarks are coupled to a temporal, spatially con- 
stant background gauge field Aq, which is represented 



in terms of the Polyakov loops, via a covariant derivative 
JP = (j) — i7oA) . A gauge coupling has been absorbed in 
the gauge fields. 

The second part in the Lagrangian, C rn repre- 
sents the purely mesonic contribution with the field 
4> = T a (a a + iir a ) and reads 

C m = Tr (d^d»<f>) - m 2 Tr(>t0) - \, [Tr(^)] 2 

-A 2 Tr (<^) 2 + c (det(<£) + det(0 t )) 

+ Tr[tf(0 + 0t)] . ( 2 ) 

The last term in Eq. ([2]) breaks chiral symmetry ex- 
plicitly where H = T a h a has in general nine external 
parameters h a whereof only two are nonvanishing for a 
(2 + l)-flavor symmetry breaking pattern. The U(1)a 
symmetry is broken explicitly by the 't Hooft determi- 
nant with a constant strength c. 

The last part of Eq. ([I]) represents the effective gluon 
field potential in terms of the Polyakov-loop variables <& 
and $ [30] which approximate the dynamical glue sec- 
tor of QCD. For vanishing chemical potential the bulk 
thermodynamics of lattice Yang-Mills is reproduced up 
to twice the deconfinement phase transition tempera- 
ture [221 [3TJ [32] . Several explicit choices of the Polyakov- 
loop potential are known in the literature [28 , 33 , 34 , see 
|29j for a comparison. In this work we employ the loga- 
rithmic version 1331 
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a{T) = a + ai 
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and the parameters ao = 3.51, a\ = —2.47, a 2 = 15.2 and 
63 = —1.75 29J. Originally the parameter Tq is fitted to 
the critical temperature To = 270 MeV of the first-order 
deconfinement transition of the quenched SU(3) lattice 
gauge theory. 

However, in full dynamical QCD the effective 
Polyakov-loop potential has to be replaced by the cor- 
responding QCD Yang-Mills (YM) part. Within a func- 
tional renormalization group approach the critical tem- 
perature of the quenched YM system could be deter- 
mined recently [35] . Because of the matter fluctuations 
the glue propagation is further dressed (see also [H])- A 
first phenomenological hard thermal loop estimate for 
this flavor and chemical potential dependence of the crit- 
ical temperature has been given in [36] and within per- 
turbation theory (see [37]). These considerations demon- 
strate that low-energy models such as the PQM model 
can be understood as a specific approximation of QCD 
and can be systematically enhanced towards full dynam- 
ical QCD. 
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To be more specific: the matter backreaction to the 
gluon sector leads to Nf and \i modifications in T , i.e., 
T T (N f ,fi) as introduced in [35] (see also [8"1 l2l? ] IB"? ] 
|3"5] for more details). 

Explicitly, we employ 

T (N f ,fj,) = T^-i/^oKW,,/.)) (4) 

where T r = 1.77 GeV denotes the r scale and ao = a(A) 
the gauge coupling at some UV scale A. The /u-dependent 
running coupling reads 

b(N f ,v) = b(N f )-b^ , (5) 

with the factor b^^^Nf. 

Finally, the grand potential in mean-field approxima- 
tion (MFA) is the sum of three terms consisting of the 
meson, U, the fermion-loop contribution Q qq , evaluated 
in the presence of the Polyakov loop and the pure gauge 



field contribution, the Polyakov-loop potential U: 

n = U(a x ,a y ) + rtg q (a Xl a y ,^^) +U ($,$) . (6) 

$ and $ denote the real Polyakov-loop expectation values 
and <j x and a y are the light and strange chiral conden- 
sates in the rotated nonstrange-strange basis (see [35] for 
further details). The temperature and density behavior 
of the four order parameters is obtained by minimization 
of the grand potential Eq. Explicitly, the mesonic 
part reads 

ui 2 c 
U(a x ,a y ) = — {a 2 x + a%) - h x a x - h y a y - ^-^ a l a y 

with six parameters m 2 ,c,\i,\2,h x ,h y and the quark- 
antiquark contribution in the presence of the Polyakov 
loop accordingly 



f=u,d,s J \ ' 



N c E qJ +T\n 



1 + 3($ + $ e - ( - E '>-f~ f *f )/T )e-( E -''f- fl f /T + e - 3{E "-f-^f )/T 
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1 + 3($ + <S> e - {E o-f + ^f )/T )e~ {E o-f + ^f )/T + e - 3{E o-f +l *f )/T | (8) 



with the flavor-dependent single-particle energies 

EqJ = y]j? + m) (9) 

and quark masses mi = ha x /2 and m s = ha y /y/2. 
The first term in quark/antiquark-loop contribution Q qq , 
Eq. ([8]), represents the ultraviolet divergent vacuum con- 
tribution. It is neglected in the no-sea MFA but can be 
renormalized with, e.g., the dimensional regularization 
scheme yielding 

= E -Mx) (10) 

/— u,d,s 

where the regularization scale A has been introduced. 
Note that the A dependence is completely absorbed by 
the mesonic potential U. Hence all thermodynamic quan- 
tities are finally independent of the choice of the regu- 
larization scale at the mean-field level. With this term 
fermion fluctuations are taken into account whereas me- 
son fluctuations are still ignored. A careful investigation 
of the influence of this vacuum term in two-flavor model 
studies can be found in [4"0H4"2"] and recently for a (2+1)- 
flavor PQM model in [29] . To disentangle physical effects 
driven by fluctuations we will compare no-sea mean-field 
results with the renormalized ones where this term has 
been taken into account. 



The parameters of the model, Eq. Q, are chosen 
to reproduce experimentally known quantities like the 
pion and kaon decay constants and the scalar and pseu- 
doscalar meson masses. The only insecure and not pre- 
cisely known quantity is the sigma mass which influ- 
ences the phase structure [IS]- To compare with pre- 
vious works [39] we choose m CT = 600 MeV for the no-sea 
mean-field approximation and m a = 400 MeV for the 
renormalized model (see also the next section). 

Neglecting the Polyakov-loop potential in Eq. ^ and 
the background gauge field A in the covariant derivative 
which corresponds to $ = $ = 1 in the quark contri- 
bution Eq. ^ results in the corresponding three-flavor 
quark-meson model. Since the Polyakov loop decouples 
from the QM sector in the vacuum no parameter adjust- 
ments are necessary. 



III. PHASE STRUCTURE AND THE CRITICAL 
ENDPOINT 

In Fig. [TJthc (2 + l)-flavor chiral phase diagrams for 
the QM and the PQM models with and without the mat- 
ter backreaction via Tq([i) are presented in no-sea mean- 
field approximation [1(a)] and in the renormalized models 
[1(b)]. The crossover in the strange sector which takes 
place at higher temperatures is not shown. The almost 
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FIG. 1. Three-flavor chiral phase diagrams [upper lines, PQM models for a constant To (top line) and for To(^t); lower lines, 
QM models], (a) No-sea MFA with m a = 600 MeV; (b) renormalized models with m a = 400 MeV. The peak in the temperature 
derivative of the Polyakov loops $, $ is also shown (light dashed lines). 



degenerated positions of the peaks of the temperature 
derivative of both Polyakov loops $ and <E> are indicated 
by the light dashed lines. 

In all models the vacuum term in the grand potential 
moves the location of the endpoint away from the T axis 
towards the (i axis whereas the Polyakov loop increases 
the chiral transition temperature at [i = and hence 
moves the CEP back to higher temperatures. For ex- 
ample the PQM model in no-sea approximation without 
the matter backreaction for a constant To = 270 MeV 
yields a CEP at (T c , /z c ) = (186, 167) MeV which changes 
to (T c , /x c ) = (172, 169) MeV if the matter backreaction 
is taken into account. Without the Polyakov loop the 
CEP lies at (T c , fi c ) = (92, 221) MeV. With vacuum fluc- 
tuations the location of the CEP changes significantly to 
(T c ,Hc) = (90 (83), 283 (280)) MeV in the renormalized 
PQM model without (with) matter backreaction. In the 
renormalized QM model the CEP lies close to the jit axis 
at (T c ,fi c ) = (32, 286) MeV. 

The location and hence the existence of a CEP in 
the phase diagram is sensitive to the sigma meson mass 
m a [jUl E3] • A larger m CT in general pushes the CEP to- 
wards the \x axis. Already for m a — 600 MeV the CEP 
disappears in the renormalized models (see also [2"9"ll39| ). 
We have adapted m a = 400 MeV in the renormalized 
models such that a common value for the chiral transi- 
tion of T x sa 205 MeV at (i = is secured as in the no-sea 
mean-field approximation. 

In the PQM models the coincidence of the chiral and 
deconfmement transitions extends to larger chemical po- 
tentials if the backreaction to the YM sector is imple- 
mented via Eq. Q. Thus, the inclusion of fluctuations 
and the matter back-coupling to the YM sector reduces 
a possible chirally symmetric and confining region in the 
phase diagram as found in a large- N c limit study [44] . A 
similar result is obtained if one goes beyond mean field 
with, e.g., functional techniques such as the functional 
renormalization group [5] or Dyson-Schwinger equations 



[TO] where both transitions coincide over the whole phase 
diagram. Hence, the inclusion of the \i dependence in T 
in the Polyakov-loop potential represents a further step 
beyond the standard mean-field approximation. How- 
ever, at vanishing temperature all transitions collapse to 
the same critical \i and the Polyakov loop does not in- 
fluence the chiral sector. Of course, in this domain other 
degrees of freedom like baryons, which are usually ne- 
glected in such (P)QM- or (P)NJL-type model studies, 
become more important which we do not address here. 



A. Anatomy of the critical region 

The knowledge of the size of the critical region around 
the CEP in the phase diagram is important for the exper- 
imental detection of the endpoint in heavy-ion collisions. 
The region is obtained as a projection of different con- 
stant ratios 



of the quark-number susceptibility \q normalized with 
the free susceptibility Xq ee on to the (T, ^)-plane near 
the CEP. In Fig. [2] three ratios R, obtained in no-sea 
MFA, are plotted as contours relative to the endpoint. 
Compared to the QM model (right panel) the Polyakov 
loop compresses the critical region, in particular in the 
T-direction. With the matter back-coupling (solid lines) 
this effect is weaker and the region is twice as large in the 
T-direction. In /i-dircction the modification is even less 
pronounced. Without the matter back-coupling the chi- 
ral and deconfmement transition deviate already in the 
vicinity of the CEP. As a consequence the quark deter- 
minant is more suppressed near the chiral transition and 
hence the size of the critical region. 

The effect of the vacuum contribution in the grand po- 
tential on the size of the critical region is shown in Fig. [3] 
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FIG. 2. The critical region around the CEP for three different three-flavor models in no-sea MFA: PQM model with a 
logarithmic Polyakov-loop potential and constant and running To (a) and the QM model (b). See text for details. 
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FIG. 3. The critical region similar to Fig. [2] for the corresponding renormalized models. 



The incorporated fluctuations increase the critical region 
in direction perpendicular to the crossover line in the 
phase diagram. This effect is more prominent in the QM 
model without the Polyakov loop [3(b)]. The Polyakov 
loop affects the critical region less than in the no-sea ap- 
proximation. However, the region in the renormalized 
PQM models [3(a)] is larger than in the renormalized 
QM model. This can be understood by the observation 
that the Polyakov loop modifies the quark determinant 
mostly at moderate chemical potentials. In the renor- 
malized models the CEP is located at larger chemical 
potentials around /i ~ 280 — 285 MeV. Thus, the en- 
hanced critical region with the Polyakov loop results also 
from the higher critical temperature [cf. also Fig. 1(b)] . 
The matter backreaction modifies the critical region in a 
similar fashion as in the no-sea approximation. 

In general, fluctuations wash out the phase transition 
and the chiral crossover is much smoother. As a conse- 
quence, the size of the critical region becomes broader in 
direction perpendicular to the extended first-order tran- 
sition line. But due to the shift of the endpoint away from 
the T-axis the size of the critical region turns out to be 



smaller along the phase boundary compared to the no- 
sea MFA. Furthermore, for larger chemical potentials and 
smaller temperatures the influence of the Polyakov loop 
becomes insignificant and hence the size of the critical re- 
gion becomes comparable in both renormalized models. 

All critical regions exhibit an elongation along the di- 
rection parallel to the first-order phase transition line. 
The peak of the quark-number susceptibility follows the 
crossover line and diverges at the CEP. The divergence 
is described by a power law within the critical region. 
The shape of the critical region is induced by the path 
dependency of the corresponding critical exponents near 
the singularity as already pointed out in [45 . One rea- 
son for the elongation of the critical region in a direction 
parallel to the first-order line is that the corresponding 
critical exponent is larger US] ■ 

This behavior is elucidated in the following: In Fig. [4] 
a generic plot of the region near the CEP in the phase 
diagram with three different paths towards the CEP, la- 
beled with A, B and C (A', B' and C) for i > (i < 0) 
where i = (T — T c ) /T c is shown. Path A (A') is cho- 
sen such that it is parallel to the (extended) first-order 
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FIG. 4. Schematic plot of three different paths (labeled with 
A, B and C and with A', B' and C) towards the CEP in 
the phase diagram. The yellow-shaded inner regions denote 
the area parallel to the (extended) first-order line (dashed) 
where the critical exponent of the quark-number susceptibility 
deviates from the one outside this area. See text for details. 

line. It corresponds to the Ising model temperaturelike 
direction with a vanishing " magnetic field" component in 
the sense of Griffiths and Wheeler 05J. The remaining 
paths B (B')and C (C) do have some additional "mag- 
netic field" contributions. For the critical exponent e, 
which is defined by 

Xq ex r~ e (12) 
with the reduced distance 

r = \J i 2 + (i 2 (13) 

from the CEP with ft, = (p— // c )//x c , two different scal- 
ing regions can be distinguished. The explicit form and 
size of these regions are not universal. They depend on 
the used thermodynamic variables and the underlying 
equation of state; however, the mapping is analytic |45j . 
When the T, /j, variables are used, the inner region is very 
small. For the outer and larger region, perpendicular to 
the extended first-order line, the mean- field value e = 2/3 
is found. Within a narrow area, parallel to this first- 
order line, a value of 7 = 1 is seen which corresponds to 
the Ising temperature-like exponent for vanishing exter- 
nal field 7 = epS = 1 > e (cf. Fig. [4]). The critical ex- 
ponent does not depend on whether one approaches the 
critical point from the side with t > or from t < 0. The 
scaling behavior is confirmed and demonstrated in Fig. [5] 
where the logarithm of the quark-number susceptibility 
for three different paths towards the CEP (for t > 0) as 
a function of the logarithm of the reduced distance r is 
presented. 

In the figure the scaling of the susceptibility for the 
PQM model without Tq(/x) corrections (left panel) and 



Path 


Jo = 270 Me V 




t > 


0.669 ±0.002 


0.669 ± 0.002 


i < 


0.662 ±0.003 


0.663 ± 0.002 


A > 


0.667 ±0.003 


0.669 ± 0.002 


A < 


0.663 ±0.002 


0.662 ±0.003 


A' 


1.00 ±0.02 


1.00 ±0.02 


A 


1.01 ±0.02 


1.01 ±0.02 



TABLE I. Critical exponents of the quark-number suscep- 
tibility in the PQM model with a logarithmic Polyakov-loop 
potential for a constant and a running To for different paths 
towards the CEP. 

for the QM model (right panel) over several orders of 
magnitude is shown. The slope corresponds to the critical 
exponent. For the path B in the QM model one nicely 
recognizes the scaling crossing from the inner region with 
an exponent 7 = 1 to the outer region with an exponent 
e = 2/3. 

A similar scaling crossing in the slopes of the corre- 
sponding PQM curves is also visible in Fig. [5] with a small 
numerical uncertainty for path B very close to the CEP. 
The distance r where the crossing of the different scaling 
regimes in the figure takes place depends on the chosen 
angle of the path relative to the CEP coordinates. This 
crossing determines the boundary of the narrow yellow- 
shaded regime in Fig. [4] Summarizing the findings for 
both models, with and without the Polyakov loop, the 
inner regime tapers with a very small angle below 0.1 
degrees. 

The criticality is also not affected by the backreaction 
of the matter fluctuations to the YM sector. In Tab. U the 
critical exponents for paths parallel to the /i and T axes 
towards the CEP with and without the T (/z) corrections 
are listed. We find for these paths the expected mean- 
field critical exponent 2/3 [47]. For the paths parallel 
to the first-order transition line (labeled with A and A') 
an exponent 7 = 1 is obtained. If one goes beyond the 
no-sea mean-field approximation using the renormalized 
potentials, the critical exponents remain mean field since 
the universality is not modified. 



IV. NON-GAUSSIAN FLUCTUATIONS 

In a realistic heavy-ion collision the detection of the 
critical region, in particular of an endpoint, is hampered 
by a finite correlation length as argued in the introduc- 
tion. When we approach a critical point more and more 
non-Gaussian components in the probability distribution 
of the corresponding order parameter are developed. The 
corresponding higher moments are much more sensitive 
on a diverging correlation length and thus seem to be a 
more appropriate tool to guide the experiment to locate 
an endpoint [48] , 

Fluctuations of conserved charges are quantified by cu- 
mulants in statistics which are directly related to gener- 
alized susceptibilities [2J |4T)H51"] . They can be defined 
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FIG. 5. The scaling of the quark-number susceptibility \q- Three different paths (cf. Fig. [4| towards the CEP are shown for 
the PQM model (left panel) and for the QM model (right panel) in no-sea MFA. 



as derivatives of the logarithm of the partition function 
with respect to the appropriate chemical potentials. For 
three different chemical potentials we have accordingly 

I QUi Qn 3 gn k 

Xm, nj ,n k = yfd(m/T) n ' d{^/T)^ d(fi k /T) n " l ° gZ ' 

(14) 

where rij, . . . denotes the order of the derivatives and the 
indices (i,j, k) — (u, d, s) the quark flavor. 

The generalized susceptibilities evaluated at vanishing 
Hf are the Taylor expansion coefficients of the pressure 
series in powers of u f/T. The coefficients and the conver- 
gence properties of the Taylor series were recently ana- 
lyzed in a similar model study [52] [53]. Moreover, higher- 
order susceptibilities are also known up to eighth order 
from QCD lattice simulations [oTl I55| . 

For the corresponding nth order moment of the quark- 
number fluctuations for one uniform quark chemical po- 
tential, we write Xn- The ratio of different order suscep- 
tibilities is volume independent and can thus be directly 
linked to moment products of various baryon number dis- 
tributions. Examples are the net-proton kurtosis re and 
the variance a which can be combined to reer 2 = Xi/i^Xz) 
and the net-proton skewness S which can be expressed 
as Scr — X3/(3x2)- The factors arise because Xn are 
the quark-number susceptibilities. These moments es- 
tablish a deeper connection between theoretical predic- 
tions and experimental measurements: only in the vicin- 
ity of the CEP might the products reer 2 and Sa show 
large deviations from its Poisson value. In the crossover 
regime at smaller chemical potentials, they are close to 
1. Presently, these quantities are the subject of ongo- 
ing discussions [32] |55] [57J . Unfortunately, the existing 
statistics on the experimental side are not yet sufficient 
to draw any reliable conclusions which might point to the 
existence of a critical point in the phase diagram [58] . 

In the following we denote the nth to mth order mo- 
ment ratio as 

Rn.rn = Xn{T, H)/Xm{T, fi) , (15) 

such that, e.g., na/S — -R4.3/3. The direct model evalu- 



ation of R n .m becomes quickly cumbersome and tedious 
already for lower orders because the order parameters 
depend implicitly on the chemical potential. 

The evaluation of the ratios can be automated by 
means of a novel numerical derivative technique, based 
on the algorithmic differentiation (AD) idea. With this 
technique derivatives of arbitrary order at machine pre- 
cision can be obtained. Compared to other differentia- 
tion techniques such as the standard divided differentia- 
tion method or symbolic differentiation, the AD is faster 
and produces truncation-error-free derivatives of a given 
function. Details and a comprehensive introduction to 
the AD technique is given in [55]. Physics applications 
of this technique can be found at various places: The AD 
generalization to implicitly dependent functions which is 
also needed in this work was done in [25] and in [32] the 
Taylor coefficients to very high orders for a three-flavor 
model have been calculated for the first time at vanishing 
chemical potential. These findings confirm impressively 
the power and usability of this technique. In this work we 
extend this method to nonvanishing chemical potentials. 

We begin with a discussion of the ratios R n .2 with even 
n > 4. Since the moment X2, i-e., the quark- number sus- 
ceptibility, is always positive and continuous in the chiral 
crossover regime it cannot modify the sign structure of 
R n< 2 near the phase transition. Thus, the ratios R n 2 are 
a suitable quantity to inspect the change of sign of all 
higher-order moments Xn in the phase diagram. 

The ratio R4 2 measures basically the quark content of 
particles carrying baryon number. For vanishing quark 
chemical potential this ratio (if quark-number suscepti- 
bilities are used) tends to nine for low temperature and 
for high temperature to 6/ir 2 which follows immediately 
from the Stefan-Boltzmann limit. In a hadron resonance 
gas (HRG) model the ratio is temperature independent 
and all moments stay positive since the HRG has no 
phase transition with singularities. Hence, higher-order 
moments can differ significantly from a HRG calculation 
along the chemical freeze-out line even if the lower-order 
moments agree with HRG results and the thermodynam- 
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FIG. 6. The ratio i?4,2 in the PQM model with a constant 
T . 



ics and the particle yields at low temperature are well 
described by the HRG model [221 S31 ED] • Any deviation 
from the HRG model result might be an indicator for a 
real critical phenomenon (in equilibrium) and thus serves 
as a theoretical baseline for the analysis of heavy-ion col- 
lisions [221152] . 

In Fig. [6] the ratio i?4 t 2 calculated in the PQM model 
in no-sea MFA with a constant To around the crossover 
region is shown. The contour lines projected onto the 
temperature and chemical potential plane are close to the 
chiral crossover line in the phase diagram. One clearly 
recognizes the development of a negative region for non- 
vanishing chemical potentials. These three-flavor results 
differ from a two-flavor mean-field calculation with and 
without the vacuum term [41] . Without the vacuum term 
the two-flavor kurtosis develops a peak at [i — near the 
crossover which is a remnant of the first-order transi- 
tion in the chiral limit [61 . The fermion vacuum term 
changes the transition to second-order in the chiral limit 
consistent with universality arguments [62, 63J. For three 
massless flavor a first-order transition always emerges in 
the chiral limit with or without the vacuum term. The 
peak of the kurtosis for physical pion masses is already 
less pronounced in the no-sea mean-field approximation. 
This suppression at fi = is caused by the strange quarks 
and not by fluctuations (cf. Fig. 7 in 29 ) . Including the 
vacuum term the peak in the kurtosis finally vanishes 
also for the three-flavor case. 

The higher moments start to oscillate within a nar- 
row temperature interval for temperatures near the chiral 
transition |32j . The structures of all moments at fi = 
are related to each other and the behavior including the 
amplitudes of \n can be estimated by the temperature 
derivative of Xn-2 [53 . In contrast to the HRG model 
the ratios in our calculation become negative in the vicin- 
ity of the crossover for nonvanishing chemical potential 
as shown in Fig. [6] Since the fourth-order moment Xi 



is unaffected by chiral critical phenomena at vanishing 
chemical potential it remains positive for /x < 20 MeV. 
This changes for larger chemical potentials since then 
nontrivial contributions from higher-order moments eval- 
uated at [i = emerge and induce the change of sign in 
the ratios [5TJ|64]. This can be seen explicitly by a Taylor 
expansion of the corresponding moment around fi = 0. 

Higher-order ratios behave in a similar way [32|, 15 lj . In 
Fig.[7]the negative regions of several ratios R n , 2 along the 
chiral crossover line are compared to each other for two 
models in no-sea MFA. All negative regions are closely 
correlated with the crossover curve. For ratios with n > 4 
the negative regions are shifted slightly away from the 
crossover curve in the hadronic phase and all regions con- 
verge exactly at the CEP. The Polyakov loop sharpens 
the size of the negative region in the phase diagram. 

The effect of the vacuum fluctuations on the negative 
regions is demonstrated in Fig. [8] where similar to the 
previous Fig. [7] the results of various even R n ,2 ratios 
with the corresponding renormalized models are shown. 
The CEP is pushed towards higher chemical potentials 
and the crossover is washed out by fluctuations. The 
curvature of the crossover line seems not to be changed. 
The effect of the fluctuations is drastic in the renormal- 
ized quark-meson model [8(b)] but not so large when the 
Polyakov loop is considered [8(a)]. 

We conclude that the negative fluctuations can surely 
be attributed to critical dynamics. Furthermore, these 
findings underline once more the importance of fluctua- 
tions. Generically, one observes that all regions calcu- 
lated in the renormalized models are shifted more in the 
hadronic phase. 

The appearance of negative values in \n at fi — has 
already been suggested as a criterion to determine the 
chiral critical temperature with Taylor expansion meth- 
ods in lattice simulations. For example, the y, = values 
of Fig. 7(a) correspond to the ones in Fig. 5 of Ref. |53) . 



However, the critical temperature estimated with this cri- 
terion deviates strongly from the chiral critical one. The 
knowledge how the negative regions evolve towards the 
endpoint might be used to construct new criteria to im- 
prove the T c estimate from a Taylor expansion around 
fi = 0. 

We close this section with a discussion of R n .2 for odd 
integers n which have not been addressed in the litera- 
ture so far. They vanish at /i = for all temperatures due 
to the CP symmetry of the QCD partition function and 
change sign at the endpoint. The quark- number density 
n q /T 3 — xi is always positive for nonvanishing chemical 
potential. Hence, negative values can only be found at fi- 
nite fx/T for n > 3. The emergence of the first zero in the 
ratios is similarly related to the chiral phase boundary as 
previously discussed for the even ratios. It is instructive 
to define the distance AT = T n — T x of the first zero in 
temperature direction of the ratio R n ,2 to the crossover 
temperature T x . In Fig. M the distance AT of four differ- 
ent ratios R n ^ is plottea as a function of ji/T in no-sea 
MFA. All curves in the figure start at /i/T « except the 
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FIG. 7. Negative regions of the ratios R n ,2 along the crossover line in the phase diagram. 
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FIG. 8. Negative regions similar to Fig. [7]for the corresponding renormalized models. 



one for the ratio R^,2 which is positive at small chemical 
potentials. All curves obey a minimum whose location 
and depth depend on the model and used approximation. 
This is clearly visible in the next Fig. [10] where the vac- 
uum term has been included. For completeness we have 
also added some even ratios in the figure which behave 
very similar to the odd ratios. Thus our considerations 
are not restricted to odd ratios. 

For the PQM model [10(a)] the negative regions are 
closer to the crossover line than for the QM model [10(b)]. 
All curves converge to the CEP. All distances are nega- 
tive for n > 5 meaning that the first zero in the ratio 
Rn,2 is pushed in the hadronic phase. The ratio -R32 
stays close to the crossover line and becomes negative 
only in the vicinity of the CEP, in particular when the 
vacuum term is included, cf. Fig. [lOj For this case the 
minimum is deeper by almost a factor of 10 due to the 
smoothening of the crossover driven by fluctuations. On 
the other hand, the Polyakov loop sharpens the transi- 
tions and consequently the minimum is not as deep (a 
factor 4 smaller) as in the QM model. 

An interesting observation is the almost linear behavior 
of AT for intermediate li/T ratios in both models with 



and without the vacuum term and for all ratios. The lin- 
ear extrapolation of AT from intermediate fi/T ratios to 
larger ratios where AT vanishes might serve as an esti- 
mator for the proximity of the thermal freeze-out to the 
crossover line and the existence of a possible endpoint 
can be ruled out for smaller Li c /T c ratios. Of course, due 
to some nonlinearities near the CEP this is only a lower 
bound for fi/T. The deviations from the linear behav- 
ior are yet smaller in the more realistic PQM model in 
contrast to the QM model (cf. Fig. 10). This estimate 
could be strengthened by considering only the difference 
of the subsequent roots in the odd ratios which is inde- 
pendent of the chiral crossover temperature. In Fig. [TT] 
the relative temperature distance At = T n+ 2~T n for sev- 
eral odd integers is shown as a function of \ijT without 
the vacuum term. The curves exhibit a similar behavior 
as in Fig. [9] except for the ordering of the curves with 
respect to n. With increasing n the distance between 
subsequent ratios decreases signaling a possible conver- 
gence of the negative regions in the phase diagram. An 
estimation of the lower bound of the CEP with a linear 
fit between fj,/T = 0.4 . . . 0.7 (= 1.0 .. . 1.5) for the PQM 
(QM) model in no-sea approximation yields already a 
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(a) Polyakov-quark-meson model with running To(fi) (b) Quark-meson model 

FIG. 9. Distance AT = T n — T x of the first zero for various ratios R n ,2 to the chiral temperature T x as a function of (x/T in 
no-sea MFA. 




critical [jl c /T c k, 0.9 (2.0) for the lowest ratios (n = 3) 
which is quite close to the actual values /x c /T c 1.0 (2.4). 
For higher n these estimates become even better. Trans- 
lating these numbers to the critical temperature yields 
T c w 175 MeV (100 MeV) already for the lowest ratio 

-^3,2- 



V. SUMMARY 

In this work we addressed the role of fluctuations in 
three-flavor chiral (Polyakov-)quark-meson models and 
discussed their phenomenological consequences. The re- 
sults obtained in no-sea mean-field approximations where 
an ultraviolet divergent part of the fermion-loop contri- 
bution to the grand potential is neglected were compared 
with findings of the corresponding cutoff-independent 
renormalized models. The influence of the vacuum term 
on the critical endpoint in the phase diagram, the size 
of the critical region around the endpoint and its crit- 
icality with respect to its effect on the universality in 
general were elucidated in detail. Because of the path 



dependence of the critical exponent towards the critical 
endpoint, the critical region which we define as a projec- 
tion of constant normalized quark-number susceptibili- 
ties on the (T, /x) plane is elongated in the direction of 
the crossover line. In the renormalized models a broad- 
ening perpendicular to the crossover line and simultane- 
ously a shrinking in the direction of the chemical poten- 
tial axis were found. The broadening is reduced if the 
Polyakov loop is taken into account. The Polyakov loop 
also shifts the crossover at vanishing chemical potential 
and consequently the critical endpoint to higher temper- 
atures. On the other hand, in the renormalized models 
the location of the endpoint is pushed away from the 
temperature- towards the chemical potential axis which 
is a genuine effect of fluctuations. This shift can be com- 
pensated by reducing the experimentally insecure sigma 
meson mass in these models accordingly. In addition, 
the matter backreaction to the pure Yang-Mills sector in 
the Polyakov-quark-meson models which represents a fur- 
ther step beyond the mean-field approximation has also 
been implemented and the impact on the location of the 
endpoint and its critical region was analyzed. With the 
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matter back-coupling improvement the chiral light phase 
transition and the peak in the temperature derivative of 
the Polyakov loops coincide over a larger region in the 
phase diagram when we synchronize both transitions at 
vanishing chemical potential by fixing the Polyakov-loop 
parameter T . This coincidence increases even further 
when more fluctuations are taken into account as seen 
in recent functional renormalization group and Dyson- 
Schwinger calculations. 

The scaling properties of the quark-number suscepti- 
bility in the vicinity of the endpoint are not modified 
by the fermion-loop vacuum fluctuations and a crossover 
phenomenon between two different scaling regimes was 
clearly seen. These regimes also determine the explicit 
shape of the critical region in the phase diagram which 
is not universal property but depends on the used ther- 
modynamic variables and underlying equation of state. 

As argued in the introduction, higher-order cumulants 
or generalized susceptibilities are more sensitive on the 
diverging correlation length and might be an appropri- 
ate quantity for the experimental search of an endpoint 
in the phase diagram. They are given as corresponding 
chemical potential derivatives of the partition function 
and are thus related to the Taylor expansion coefficients 
of the pressure series. 

By applying a novel numerical derivative technique, we 
could calculate higher moments to very high orders at fi- 
nite chemical potential. In contrast to the HRG model 
which does not have a phase transition with singulari- 
ties, the kurtosis stays positive. In our models and in- 
dependent of the used approximation the higher-order 
moments differ significantly from the HRG model result 
along the freeze-out curve due to the existence of an end- 
point. As a consequence the ratios between various mo- 
ments develop a negative region in the phase diagram 
close to the chiral transition line. Fluctuations smoothcn 
the transitions and hence the ratios and the negative re- 
gions are broader and are shifted towards the hadronic 
phase in the renormalized models. All regions converge 



at the endpoint. 

In order to quantify this general trend we have in- 
troduced two new quantities: the distance of the first 
zero in temperature direction of various ratios to the 
crossover temperature and the relative distance between 
subsequent roots in the ratios which does not require the 
knowledge of the insecure chiral crossover temperature. 
Both quantities have a minimum as a function of the 
fi/T ratio. By using a linear extrapolation for interme- 
diate /i/T ratios as an estimator we could rule out the 
existence of an endpoint for smaller [ijT ratios. This 
estimator might be useful for lattice simulations where 
only the negative regions for the first few ratios for small 
densities might be available. 

An interesting difference to a corresponding two-flavor 
PQM calculation concerns the ratio at vanishing 

chemical potential: without the vacuum term the two- 
flavor kurtosis develops a peak at /x = near the 
crossover. The vacuum term smoothens the transition 
and the peak disappears finally. In the three-flavor cal- 
culation the peak is already suppressed in the no-sea ap- 
proximation and hence is not driven by fluctuations. 

It would be interesting to investigate higher moments 
of the so far neglected electric charge or strangeness fluc- 
tuations in particular beyond the mean-field level with 
the functional renormalization group. Work in this di- 
rection is ongoing. 
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